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Abstract 

, We study effects of super symmetric particles in various rare B decay processes as well as in 

' the unitarity triangle analysis. We consider three different supersymmetric models, the minimal 

. supergravity, SU(5) SUSY GUT with right-handed neutrinos, and the minimal supersymmetric 

, standard model with U(2) flavor symmetry. In the SU(5) SUSY GUT with right-handed neutrinos, 

cn 

' we consider two cases of the mass matrix of the right-handed neutrinos. We calculate direct and 

mixing-induced CP asymmetries in the 6 — > 57 decay and CP asymmetry in B^i 4>Ks as well 
Ph! as the Bs~Bs mixing amplitude for the unitarity triangle analysis in these models. We show that 



large deviations are possible for the SU(5) SUSY GUT and the U(2) model. The patterns and 
correlations of deviations from the standard model will be useful to discriminate the different SUSY 



X 

. models in future B experiments. 
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I. INTRODUCTION 



Success of B factory experiments at KEK and SLAC indicates that B physics is important 

to explore origins of the flavor mixing and the CP violation in the quark sector. The CP 

asymmetry in the Bd —>■ J/ipKs mode is precisely determined, and a CP violation parameter, 

sin 201 (or sin 2/3) is found to be consistent with the prediction in the standard model (SM) 

Q]. In addition, branching ratios and CP asymmetries in various rare B decays have been 

reported. In future, we expect much improvement in measurements of CP violation and rare 

B decays at e~^e~ B factories as well as hadron B experiments P]- 

Goals of future B physics are not only to very precisely determine the parameters of the 

n 

Cabbibo-Kobayashi-Maskawa (CKM) matrix elements |3|, but also to search for new sources 
of CP violation and flavor mixings. For instance, scalar quark (squark) mass matrices could 
be such new sources in supersymmetric models. Since the flavor structure of the squark 
mass matrices depends on the mechanism of supersymmetry (SUSY) breaking at a higher 
energy scale and interactions above the weak scale through the renormalization, future B 
physics can provide quite iniportant information on the origin of SUSY breaking. 

In our previous papers 0,1^; we studied the flavor signals in three different models namely, 
the minimal supergravity (mSUGRA), the SU(5) SUSY GUT with right-handed neutrinos, 
and the minimal supersymmetric standard model (MSSM) with U(2) flavor symmetry j^, 
Q|. These models are typical solutions of the SUSY flavor problem. They have different 
flavor structures in the squark mass matrices at the electroweak scale. Thus, they may be 
distinguished by low energy quark flavor signals. We calculated SUSY contributions to the 
Bd-Bd, Bg-Bs, and K^-K^ mixings in these models, and showed that the consistency test of 
the unitarity triangle from angle and length measurements are useful to discriminate these 
models. 

In addition to the consistency test of the unitarity triangle, there are several promising 
ways to search for new phvsics effects through B decay processes. As pointed out in the 
context of SUSY models P, 0, comparison of time-dependent CP asymmetries in Bd — >■ 
J/ipKs, Bd (pKs, and Bd Tj'Ks provides us with information on new CP phases in the 
b s transition, because these asymmetries are expected to be the same in the SM. Recent 
results on the CP asymmetry in Bd — > ())Ks by Belle and BaBar collaborations indicate a 
2.7cr deviation from the SM prediction |l1J|. This anomaly may be attributed to new physics. 
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SUSY with R parity U, 1^, SUSY without R parity or other models Q]. Moreover, 
the CP a=yo,.et™= . the , s, p.oce. QQ have been exten.ve,y =tud.ed ™, 
models, and they would exhibit substantial deviation from the SM in some models |12|]. The 
branching ratio and decay distributions of 6 — > sll have also examined by several authors 
in different contexts of new physics, and they may probe different aspects of new physics 



lia 



19|. 

In this paper, we extend our previous analysis ^ to rare B decay processes, especially 
b ^ s transitions. We investigate the direct CP asymmetry in 6 — > 37, and the mixing 
induced CP violation in Bd Ms'~f process, where Mg is a CP eigenstate hadron with 
strangeness, and CP asymmetry in Bd (pKs process. We show that the above three 
models exhibit different patterns of deviation from the SM. These observables as well as 
those studied our previous paper Q could play aurmportaut role in uew physics search 
at future B experiments such as a super B factory j20| and hadron B experiments [2]. 

The strategy of the present work and [J] is different from most of other works. For 
each model of SUSY, we calculate SUSY effects in various mixings {Bd-Bd, Bg-Bg, and 
K'^-K^) and rare decay processes, and identify possible patterns of deviations from the SM 
predictions. We then compare patterns of the new physics signals for different models. In 
this way we may be able to distinguish different models of SUSY breaking scenarios, or at 
least obtain important clues to identity the SUSY breaking sector. Most of past works deal 
with a specific observable signal in a particular SUSY model. The strategy of combining 
various information in B physics will be important in future, especially in the days of a 
super B factory and dedicated hadron B experiments. The purpose of the present work is 
to demonstrate how such global analysis in B physics is useful to explore the flavor structure 
of SUSY breaking sectors. 

In the present work, some of our calculations are re-analysis of past works under most 
up-dated phenomenological constraints, and others are new studies. Even in the case that 
a similar calculation can be found in the literature for a specific process and a model, we 
repeat the calculation in order to treat various processes in a uniform way. In the mSUGRA 
model, Bd-Bd, Bg-Bg, and K^-K^ mixing, and the direct CP violation of 6 — > 57 were 
studied previously. Although the mixing-induced CP violation in Bd Ms7 process and 
CP asymmetry in Bd —>■ <f)Ks process have not explicitly considered before, it was recognized 
that such processes did not induce interesting signals once the severe constraints from various 
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electric dipole moments were applied. We have confirmed this by explicit calculations. For 
the SU(5) SUSY GUT with right-handed neutrinos, we have found that flavor signals are 
quite different for different choices of the heavy right-handed neutrino mass matrix. The 
"degenerate case" deflned in the next section was considered in jj, 21, 22|. In particular, 
SUSY effects on the mixing amplitudes were studied in and the mixing- induced 

CP violation in Ms'j process was also considered in j2l|. The direct CP violation 

in b —>■ s'j and the CP asymmetry in <pKs process are analyzed here for the first 

time. The analysis of B physics signals in "non-degenerate case" is new. We compare the 
implications of consistency test of the unitarity triangle for degenerate and non-degenerate 
cases, which has been done only for degenerate case in |^. We also present analysis of rare 
decay processes and the Bs mixing for non-degenerate case, which partly overlaps with a 
recent work in For the U(2) model, fiavor signals has been considered previously for 
mixing amplitudes in 0, but a quantitative analysis of various rare decay modes has 
not yet appeared in the literature. Preliminary results of our analysis are presented in a 
workshop report [5], but the non-degenerate case is not included in Ref. j^. 

Although we tried to clarify a characteristic feature of each model, this cannot be com- 
plete because of large numbers of free parameters. There might be some parameter space 
where new contributions become particularly important, which are not shared in a generic 
parameter space. For example, it was pointed out recently that neutral Higgs boson ex- 
change effects may contribute to the various flavor changing process, especially for a large 
value of the ratio of the two vacuum expectation values (tan/3) [2^, but we do not take into 
account these effects. Although we estimate that such an effect is small for the parameter 
region presented numerically in this paper, this effect will be potentially important. 

This paper is organized as follows. In Sec. II, we introduce the three models. The 5^-5^ 
mixing, the Bg-Bg mixing, CP asymmetries in 6 ^ 57 and B^ (pKs are discussed in 
Sec. III. The numerical results on these observables are presented in Sec. IV. Our conclusion 
is given in Sec. V. 

II. MODELS 

In this section, we give a brief review of the models. They are well-motivated examples 
of SUSY models, and are chosen as representatives that have distinct flavor structures. A 
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detailed description of these models can be found in Ref. 



A. The minimal supergravity model 

In the mSUGRA, SUSY is spontaneously broken in the hidden sector and our MSSM 
sector is only connected to the hidden sector by the gravitation. The soft breaking terms 
are induced through the gravitational interaction, and the soft breaking terms have no new 
flavor mixing at the scale where they are induced. 

The soft breaking terms are specified at the GUT scale by the universal scalar mass 
(mo), the universal gaugino mass (M1/2), and the universal trilinear coupling (Aq). The soft 
breaking terms at the electroweak scale are determined by solving renormalization group 
equations. 

In this model, the only source of flavor mixings is the CKM matrix. New flavor mixings 
at the electroweak scale come from the CKM matrix through radiative corrections. 

As for CP violation, in addition to the CP phase in the CKM matrix, we have two new 
CP phases. One is the complex phase of the /i term (0^), and another is the phase of Aq 



(EDM) to CP 
the neutron 



{(f) a)- Since the potential sensitivity of the neutron electric dipole moment 
violations in low energy SUSY models was stressed by Ellis et al. in Ref. [2 
EDM and the electron EDM have been studied in detail by several authors in different 
ccext of low e„e.g. SUSY HflQ. TKe bo.cn U.e *e .SUGRA ,s that ^a„d 
ipA contribute to the neutron and the electron EDM's, and experimental constraints Q^Sl 
on these phases are very severe. Taking these EDM constraints into account, effects of new 



CP phases on K and B physics have turned out to be small 30 1. 



B. The SU(5) SUSY GUT with right-handed neutrinos 

In the last decade, three gauge coupling constants were determined precisely at LEP and 
other experiments, and the measured values turned out to be consistent with the prediction 
of supersymmetric grand unification. Furthermore, recent developments of neutrino exper- 
iments established the existence of small finite masses of neutrinos 

HQS, 

which can 



be naturally accommodated by the seesaw model 



3- 



Guided by these observations, we 



consider SU(5) SUSY GUT with right-handed neutrinos. 
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In this model, the soft breaking terms are the same as in the mSUGRA model at the scale 
where the soft breaking terms are induced. Unlike the mSUGRA model, the SU(5) SUSY 
GUT with right-handed neutrinos has new sources of flavor mixing in the neutrino sector. 
A large flavor mixing in the neutrino sector can affect the right-handed down type squark 



sector t 
Ref. b, 



irough GUT interactions. Quark flavor signals in this model have been studied in 
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211^,^ 



In the seesaw model, the neutrino mass matrix is written as 

{m^r = {h2)'iy.nM^Tiy.y^ , (i) 

where i/i, is the neutrino Yukawa coupling constant matrix, Mjy is the mass matrix of right- 
handed neutrinos, (/i2) denotes the vacuum expectation value of one of the Higgs fields h2, 
and i,j,k,l are generation indices. In the basis that the charged lepton Yukawa coupling 
constant matrix is diagonal, this neutrino mass matrix is related to^the observable neutrino 
mass eigenvalues and the Maki-Nakagawa-Sakata (MNS) matrix ^i^J as 

(m^r = (^MNs)^'^'(n^Ns)'^' ■ (2) 

In this model, the scalar lepton (slepton) masses and the A terms have the mSUGRA 
type structure at the Planck scale (Mp) as 

(mlr = ml5'i , {mlr^ = ml5'^ , {Ae^ = m.AoylS'^ , (3) 

where m\ and the mass squared matrices of sleptons and Ae denotes the slepton 

trilinear scalar couplings. However, at the scale M^, where the right-handed Majorana 
neutrinos are decoupled, new flavor mixings are generated by the renormalization group 
effects j^^l- III leading logarithmic approximation, they are given as 

{mlf ^ - ^mg(3 + \A,\'){yly.r In ^ , (4a) 
{mlf ~0 , (4b) 
{AeT ^ - ^,m,A,yl{yly,f In ^ , (4c) 

for i ^ j. Consequences of these mixings on lepton flavor violating processes have been 
investigated from various aspects We see that, in this model, lepton flavor violating 

processes such as — > 67 are sensitive to the off-diagonal elements of y^yu- 



We consider two cases in the SU(5) SUSY GUT with right-handed neutrinos in regard to 
the spectrum of the right-handed Majorana neutrinos. One is the case that all the masses of 
heavy Majorana neutrinos are the same (degenerate case). In this simplest case, the flavor 
mixing of the neutrino sector is only caused by Uu because there is no flavor mixing in Mat. 
Since the large mixing in the MNS matrix implies that the off-diagonal elements of is 
large, the /z — 67 branching ratio is enhanced in the wide region of the parameter space 
and exceeds the experimental bound in some parameter regions. In order to suppress the 
/i ^ 67 branching ratio, we consider an elaborated case that Mjv is not proportional to the 
unit matrix (non-degenerate case). In this case, the neutrino mixing comes from both 
and Mat. If the large mixing in the MNS matrix originates from Mat, the corresponding 
off-diagonal element of need not be large. Thus, the /z — >• 67 decay rate can be suppressed 
in this case. 

We have new CP phases in this model. They are classified in the following three classes: 

(i) The CP phases in the mSUGRA, i.e., <;/)^ and (^a- 

(ii) CP phases in the neutrino sector. There are six physical complex phases in y^ and 
Mat in the basis in which the charged lepton mass matrix is real and diagonal. From 
the combination of these six CP phases, we obtain three CP phases in the low energy 
region, i.e., one Dirac CP phase and two Majorana CP phases. 



iii) GUT CP phases |9|, |22, |40| . The quark and lepton superfields are embedded in 10 and 



5 representations of SU(5) as 



10,: 



{Q„e-^*?(0f/)i,e^^"E,} , 5, = {A,e-^<^"L,} , (5) 



where V is the CKM matrix at the GUT scale, gi(3,2, 1/6), [7:^(3,1,-2/3), 
A;(3, 1,1/3), Lj(l,2, — 1/2), and _Ej(l,l,l) are quark and lepton superfields in the 
ith generation with the SU(3) x SU(2) x U(l) gauge quantum numbers in parentheses. 
The phases 0f and 0f obey the constraints + 02 + ^3 = and 0f + 0^ + ^3 = 0. 
Before the SU(5) is broken, CP phases 0f and 0^ have physical meanings and they 
may play an important role in the flavor physics through the renormalization group 
effect above the GUT scale. 
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C. A model with U(2) flavor symmetry 



An alternative solution of the flavor problem of SUSY is introducing some flavor symme- 
try. U(2) flavor symmetry is one of such symmetries 0, Q]. We consider the model given 
in Ref. jlf. In this model, the quark and lepton supermultiplets in the first and the second 
generations transform as doublets under the U(2) flavor symmetry, and the third generation 
and the Higgs supermultiplets are singlets under the U(2). 

In order to reproduce the correct structure of the quark Yukawa coupling matrices, we 
assume the following breaking pattern of the U(2): 



U(2) U(l) l(no symmetry) . 



(6) 



With this assumption, we obtain the quark Yukawa coupling matrix yq and the squark mass 
matrices m^: 



y4 = yQ 






age' 
















1 



Q = U,D 



(7) 



m 



X 



[m. 



X\2 



\ 



,X = Q,U,D 



(8) 



1 

1 + rge^ r^e 

where e and e' are order parameters of the U(2) and U(l) symmetry breaking respectively 
and they satisfy e' <^ e ^ 1, and Yq, aq, bg, cq, dq, and r"^ are dimensionless parameters 
of 0{1). As for the squark A terms, they have the same structure as the quark Yukawa 
coupling matrices: 

/ n . n \ 



dqe' 
-aqe' dqe bqe 
cqe 1 J 



Q = U,D. 



(9) 



In general, though being of 0{1), dq, bq, cq, and dq take different values from the corre- 
sponding parameters in Eq. (0), and we expect no exact universality of the A terms in this 
model. 
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In the mass matrices of sfermions in this model, the degeneracy between masses of the 
first and the second generation is naturally realized. On the other hand, the mass of the third 
generation may be separated from the others. There exist flavor mixings of 0{e) between 
the second and the third generations of sfermions. These are new sources of flavor mixing 
besides the CKM matrix. 

There are several efforts to explain the observed neutrino masses and mixings in SUSY 



models with the U(2) flavor symmetry (or its discrete relatives) [41]. However, the purpose 
of this paper is to illustrate typical quark, especially bottom, flavor signals in several SUSY 
models and to examine the possibility to distinguish them. A detailed study of neutrino 
masses and mixings and lepton flavor signals in models with flavor symmetry is beyond the 
scope of the present work. Hence, in the following analysis, we will not consider the lepton 
sector in the U(2) model. 

III. PROCESSES 

The processes considered in the following are the 5^-5^ and the Bg-Bg mixings, b 57, 
and B (pKs- We consider the effective Lagrangian that consists of AB = 1 terms and 
AB = 2 terms: 



I — I 



^ — ^AB=l + ^AB=2 , 



(10) 



^AB=1 —C2LO2L 




(11) 
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where O's are 



O2L ={.Sc.l^CLa){.Cpl^hLp) , 

O2L ={Sal^ULo){U(3-i%Lf5) - {SaYcLa){cpYbLI3) , 



(12a) 
(12b) 



(12d) 



(12c) 



(12e) 



=\{Sa[r,Y]ha){Sf,[l„lu]sLp) , 



(12h) 




(12f) 



(12i) 



Among the above Wilson coefficients, C2L is dominated by the tree contributions from the 
SM at the weak scale and the others are induced by loop effects. Therefore, one obtains 
= tuC2Li where e„ = —V*gVub/Vfl.Vtb. The coefficients Cll, C'^lw ^'^'^ ^fk ^^^^ 
dominated by the SM contribution because of the QCD correction below the electroweak 

n 

scale. Cab=2 is described in our previous paper 

We discussed the Bd-Bd and the Bg-Bg mass splittings Attib^ and Am^^ in Ref. Q]. In 
this paper, we consider both the direct and the mixing induced CP asymmetries in b s'j 
and the CP asymmetry in 5 — ^> (pKs in addition to the above AB = 2 process. 

About the B^-B^ mixing processes, the mixing matrix elements Mi2{Bd) and Mi2{Bs) 
are defined as 



M,2{Bg) 



1 



{Bq\CAB=2\Bq) , 



(13) 



where q = d,s. We can express Am^^ and Aitlb, in terms of Mi2{Bq) as 



AruB, =2|Mi2(5g)| . 



(14) 
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The direct CP asymmetry in the inclusive decays B X^'j is defined as [1 



r(g^x,7)-r(E-.x,7) 



"3 



A\C7L\'+\Cm\ 



-Imr2lm[(l-e„)C2LC;j 



80 

—iTlm{euC2LC;L) 



+ -nlmiCsLC;^) - Im/27lm [(1 - tu)C2LC'^j\ 



+ -Im/27lm [(1 - e„)C2LC*J + (L ^ i?) 



where the functions and j^i are found in Ref. 



(15) 



42|. 



The time-dependent CP asymmetry in the decays to a CP eigenstate jcp is given by 



r(5,(t)-./cp) + r(i?,(t)^/cp) 



^cp(^rf^/cp)cosAmB,t 
+ Agi5(5rf^/cp)sinAmB,t 



(16) 



In the b ^ s'y decays, we consider the time-dependent mixing induced CP asymmetry in 
Bii Ms7, where Mg denotes a hadronic CP eigenstate which includes a strange quark such 



as K* or Ki. A^f{Bd M^'j) is given as 



3 



2Im(e-^<^^C7LC7p) 



where 



Mi2(5rf) 



+ IC 



7R 



|Mi2(5,)| 

As for (pKs, we consider 74™p(i?rf — > (pKs), which is given as 



(17) 



(18) 



(19) 



\A\^ + \A\^ 

where A and A are decay amplitudes of B^ (pK and (pK respectively. Since 

Bd —* (pKs is a hadronic decay mode, the calculation of the decay amplitude suffers from 
large theoretical uncertainties. Here we use a method based on the naive factorization. 



Details of the calculation of A is given in Refs. 
ansatz, we obtain 



y,Q. 



Using the naive factorization 



A = ~ {(j)K\C\B) 
= -Hv 



(2) 
LR 



V 



TL 



5 Ht (2) 



0,4 
47r9 



kdmCsl + (L i?) 



(20) 
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where Hy = 3{(j)K\OLL\B) and Ht = -{6/7){(j)K\0^^l\B). P^''> {q"^ , ml) comes from the 
one- loop matrix element of O-^l, is the momentum transfer of the exchanged gluon, and 
kdm is an 0{1) coefficient |43i] that parametrizes the matrix element of Os- The concrete 

brmof P(J'^(g2,m2) isP^?'\q\ml) = -{os / {47r)){G{ml, q^, ml) + 2/3) for the NDR scheme 

3 with 

G{ml,q^ml) = -4 dxx{l-x)\n{ ^ . (21) 

In the constituent quark model, Ht/ Hy is proportional to m^/ms, and thus the contribu- 
tions from C^l and C^l are neglected. 

In our calculations, we take = ml/2 and kdm = 1- Though these approximations are 
difficult to be justified in QCD, we employ them for an illustration that may provides the 
correct order of magnitude. 



IV. NUMERICAL ANALYSIS 



A. Parameters 



1. Parameters in the minimal supergravity model 

The parameters we use in our calculation are almost the same as those used in our 
previous paper P] except for CP phases. In our calculation, we treat the masses and the 
mixing matrices in the quark and lepton sectors as input parameters that determine the 
Yukawa coupling matrices. 

The CKM matrix elements Vus, Vcb, and \Vub\ are determined by experiments inde- 
pendently of new physics because they are based on tree level processes. We adopt 
Vus = 0.2196 and Vcb = 0.04 in the following calculations and vary \Vub\ within a range 
iKife/Kbl = 0.09 ± 0.01. Note that the current error of \Vub\ is estimated to be larger than 
this value but we expect theoretical and experimental improvements in near future. We vary 
the CKM phase 03 = arg(— V^*^\4d/V^'J,Vcd) within ±180°, because it is not yet constrained 
by tree level processes free from new physics contributions. 

As for the SUSY parameters, we take the convention that the unified gaugino mass M1/2 
is real. It is known that 0^ is strongly constrained by the upper bound of EDM's, while the 
corresponding constraint on 0^ is not so tight j27|. Accordingly, we fix 0^ as 0° or 180° at 

12 



the electroweak scale. We vary the universal scalar mass mo, M1/2, and the proportional 
constant of A terms to Yukawa coupling matrix AqUIq within the ranges < rriQ < 3TeV, 
< Ml/2 < ITeV, \Ao\ < 5, and -180° < (pA < 180°. We take the ratio of two VEV's 
tan/5 = {h2)/{hi) = 30 or 5. 



2. Parameters in the SU(5) SUSY GUT with right-handed neutrinos 

In the SU(5) SUSY GUT with right-handed neutrinos, we need to specify the parameters 
in the neutrino sector in addition to the quark Yukawa coupling constants given in the 
previous discussion. We take the neutrino masses as m^^ — m^^ = 3.5 x 10~^eV^, mL — — 
6.9 X 10~^eV^, and rrij^^ ~ O.OOleV, and the MNS mixing matrix as 



V2 



( 



MNS 



(22) 



CqCis SqCis Si3 

"■SoCatm ~ C0SatmSl3 C0Catm ~ ^QSatmSlS SatmCl3 
y ■SQSatm ~ CoCatm-SlS ~C0Satm ~ ■S0Catm'5l3 CatmClS j 

(cj = COS 6*4, Sj = sin6'j) with sin^ 2^atm = 1, tan^ 6*© = 0.420, and sin^ 26'i3 = 0. These mass 
squared differences and mixing angles are consistent with the observed solar and atmospheric 



neutrino oscillations [31|, the K2K experiment [3^, and the KamLAND experiment 
Only the upper bound of sin^ 26*13 is obtained by reactor experiments 4^ , and we take the 
above value as an illustration. We do not introduce the Dirac and Majorana CP phases in 
the neutrino sector for simplicity. 

We consider = 4.0 x lO^^GeV and Mr = 4.0 x lO^^GeV for the degenerate case. 
(In this case, Mn is the same as the common mass of the right-handed neutrinos.) In the 
non-degenerate case, we take the neutrino Yukawa coupling matrix as 



'^Fii ^ 



ylvu 



Y22 



\ 



(23) 



Y*, Yss J 



in order to avoid too large SUSY contribution in the branching ratio of /i — > 67. The masses 
and mixings of the right-handed neutrinos are determined to reproduce the observed neutrino 
masses and Vmns- In Table HI we show the numerical values of the neutrino parameters in the 
non-degenerate case. We integrate out the right-handed neutrinos at Mr = 4.0 x lO^^GeV 
in the non-degenerate case. 



13 



tan P 



eigenvalues of Mat (xlO^^GeV) 



0.13 
0.099 -0.099 
0.46 0.46 



4.4 
0.56 
1.7 



30 



0.13 
0.098 -0.098 
0.47 0.47 



4.5 
0.57 
1.8 



TABLE I: The neutrino parameters used in the non-degenerate case. We show the neutrino Yukawa 
couphng matrices and the mass eigenvalues of the right-handed neutrinos for tan /3 = 5 and tan (3 = 
30. These Yukawa coupling matrices give the structure as given in Eq. 

For the GUT phases 0f and 0f , we take 0f = and vary 0^ within -180° < 0f < 180° 
while 0f + 02 + 03 = is satisfied. 

The soft SUSY breaking parameters in this model are assumed to be universal at the 
Planck scale, and the running effect between the Planck scale and the GUT scale is taken 
into account. We scan the same ranges for mo, M1/2, 0a, and \Aq\ as those in the mSUGRA 
case. 

In Fig. we show the branching ratio of /i — ^> 67 in both the degenerate and the non- 
degenerate cases of the SU(5) SUSY GUT with right-handed neutrinos. As seen in this 
figure, the experimental constraint on the parameter space is more strict for the degenerate 
case than the non-degenerate case. Both in the degenerate and the non-degenerate cases, 
the SUSY contribution to ji ^ e-^ becomes larger for larger Mr 0, 0]. However, the 
contribution is less significant in the non-degenerate case for a similar Mr 2^. In fact, 
for Mr = 4.0 X lO^^GeV in the degenerate case, the SUSY contribution to /i — 67 is so 
large that most of the parameter region is excluded when tan/5 = 30. While, in the non- 
degenerate case, a large part of the parameter region is allowed. In the following, we take 
Mr = 4.0 X lO^^GeV for the degenerate case, and Mr = 4.0 x lO^^GeV for the non-degenerate 
case. 
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m(v) [ GeV ] m(v) [ GeV ] 



FIG. 1: The branching ratio of /i ^ 67 as functions of the lightest sneutrino mass in the SU(5) 
SUSY GUT with right-handed neutrinos. The dotted hue shows the experimental upper bound. 

3. Parameters in the U(2) model 

In the U(2) model, the symmetry breaking parameters e and e' are taken to be e = 
0.04 and e' = 0.008, and the other parameters in the quark Yukawa coupling matrices 
are determined so that the CKM matrix and the quark masses given in Sec. IIV A II are 
reproduced. The detailed discussion to determine the quark Yukawa coupling matrices is 
given in Ref. 0|. 
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There are many free parameters in the SUSY breaking sector as shown in Eq. (jH}. In 
order to reduce the number of free parameters in numerical calculations, we assume that 



02 U2 _D2 _ 2 



(24) 



Q U D _ 



(^j ) = (22) , (23) , (33) . (25) 

We scan the ranges for these parameters as < mo < 3TeV, — 1 < < +1, < rsa < 4, 
|r23| < 4, and —180° < argr23 < 180°. We assume that the boundary conditions for the A 
parameters are the same as the mSUGRA case for simplicity. 

B. Experimental constraints 

In order to obtain allowed regions of the parameter space, we consider the following 
experimental results: 

. Lowe. IMts on the of SUSY particles and the Higg= bosons given by direct 

searches in collider experiments |4fil |. 

• Branching ratio of the 6 — > 57 decay: 2 x lO^'' < B(6 57) < 4.5 x lO"'' 0. 

• Upper bound of the branching ratio of the /i ^ 67 decay for the SUSY GUT cases: 
B(/i ^ 67) < 1.2 X 10"" [48]. 

• Upper bounds of EDM's of the neutron and the electron: \dn\ < 6.3 x lO^^^e ■ cm 
and \de\ < 4.0 x 10~^''e ■ cm 

• The CP violation parameter ek in the K^-K^ mixing and the Ba-Ba mixing parameter 

n n 

Am^^ [49]. As for the B^-B^ mixing parameter Am^^, we take Am^^ > 13.1ps ^ |50|. 

• CP asymmetry in the B ^ J / ipKs decay and related modes observed at the B factory 
experiments [l|. 

C. Numerical results 

1. Unitarity triangle analysis 

As in our previous work j^, we search possible values of A^plB^ J/ipKs), ^ttib^, 
ArriB^, and 03 under the constraints stated above. The results for the mSUGRA and the 
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FIG. 2: AmB^/^niBj^ versus the mixing-induced CP asymmetry of Bfi J/tpKs and 03 in 
the SU(5) SUSY GUT with right-handed neutrinos. The hght-colored regions show the allowed 
region in the SM. The curves show the SM values with iKfe/Ktl = 0.08, 0.09 and 0.10. This plot 
corresponds to Fig. 5 of Ref. P|. 

U(2) model are not shown here, since they are similar as Fig. 5 in Ref. [J] apart from slight 
changes in some input parameters. 

In Fig. 121 we show the above quantities in the non-degenerate case of the SU(5) SUSY 
GUT with right-handed neutrinos for tan/? = 30 as well as the degenerate case. In the 
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mSUGRA 


SU(5) SUSY GUT 


U(2) 






degenerate 


non-degenerate 






- 


V 




V 


AmBjAniBa 




V 


V 


V 



TABLE II: Possible deviations of 03 and AniB^/AmB^ from values expected in the SM. -y/ means 
large deviation. 

degenerate case, as we found in our previous paper j^, the SUSY contribution to the 7^°- 
K'^ mixing is large and Ek is significantly affected. This is because of the large 1-2 mixing 
in the squark sector. However, the SUSY contributions to the Bd-Bd mixing and the Bg-Bg 
mixing are not important. Thus the correlation among AtubJ AniB^-, A^^p{Bd J/ipKs), 
and 03 is very similar as the SM. On the other hand, in the non-degenerate case, the 2- 
3 mixing in the squark sector is enhanced, and the 1-2 mixing and the 1-3 mixing are 
suppressed. This means that the correlation among ArriBj ArriB^, AQp{Bd — * J/ipKs), 
and 03 may differ from the SM, because of non-negligible SUSY contributions to the Bg-Bg 
mixing. 

From this figure, we see that the allowed range of 03 depends on the value of Ams^/Amg^ 
in each case. For example, if AtubJ Aitlb^ is consistent with the SM (~ 35), we observe 
03 ~ 60° in the degenerate case, and 45° ^ 03 ^ 75° in the non-degenerate case. In the 
case that AuibJ ArriB^ is larger than the SM, e.g., AniBjArriB^ ~ 55, 03 is ~ 40° in the 
degenerate case, while the allowed range is expected as 45° ^ 03 ^ 60° in the non-degenerate 
case. This indicates that a 03 measurement is important to distinguish the two cases. In 
Table HH we summarize possible deviations of 03 and Am^^ from the SM in each model. 

In the above and the following numerical calculations, as mentioned in the introduction, 
we do not take the tan /5-enhanced radiative corrections in the neutral Higgs couplings 
into account. These corrections could be significant for the B — B mixing amplitudes in 
the parameter regions of large values of tan (3 and small masses of the heavy Higgs bosons, 
because they scale as (tan (3)'^/ {mH^Y I^S]- Accordingly, it is unlikely that this effect changes 
our numerical results momentously in the region tan/3 ^30. As for B^, <pKs, it is shown 
by Kane et al. that the contribution from the neutral Higgs exchange diagrams is small. 

In Fig. El we show allowed regions in the ReMi2(-Bs) and hn.Mi2{Bg) plane for the three 
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FIG. 3: Real and imaginary parts of Mi2{Bs). The SM value is shown in a light color. 

models in the case of tan f3 = 30. In the mSUGRA, the deviation from the SM is less 
than 5%, and the SUSY contributions to the complex phase of Mi2{Bs) are negligible. In 
the non-degenerate case of the SU(5) SUSY GUT with right-handed neutrinos, the SUSY 
contribution is large in both the real and the imaginary parts of Mi2{Bs). They can be as 
large as 30% of the SM contribution to \Mi2{Bs)\. This is in contrast with the degenerate 
case. As seen in Fig. 01 the SUSY contribution to Mi2(-Bs) in the degenerate case is tiny, 
since the constraint to the parameter space imposed by the fi —>■ ej branching ratio is very 
strict. In the U(2) model, there are SUSY corrections of the order of 20% or larger to 
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TABLE III: The relation between the Wilson coefficients and the observables. ^J means that the 
coefficient gives a main contribution to the observables. 

Mi2{Bs). We have studied the case of tan/5 = 5 as well. We have found that the allowed 
regions of Mi2{Bs) are similar to those for tan (3 = 30. From the experimental point of view, 
^^Bs = 2|Afi2(-Bs)| will be measured by using Bg decays such as Bg —>■ DgTi in hadron B 
experiments j^. The phase of Mi2{Bs) can be measured by observing CP violation in Bg 
decays such as Bg J/ipcp. 

2. Rare B decays 

Here, we discuss rare B decays in the three models. We first present SUSY contributions 
to the Wilson coefficients of the dipole operators Ctl, C-jr, Csl, and Csr- In Table HTH we 
show the relation between these Wilson coefficients and observables. 

The real and imaginary parts of Cjl and Cyr at the bottom mass scale divided by the 
SM value of C^l are plotted in Fig. Hlfor tan /5 = 30. For tan (3 = 5, SUSY contributions are 
less significant, and we mainly consider the tan/? = 30 case in the following. 

In the leading order approximation, the branching ratio of 6 — 57 is proportional to 
IC'tlP + ICtkP- In the mSUGRA, however, the SUSY contributions to Cjr is very small, 
because of no new flavor violation in the right-handed squark sector. Thus, the 6 — 57 
branching ratio constrains ICjlI. In addition, the SUSY contributions to the phase of Cjl, 
which is dominated by the phase 4>a, is small due to the constraint from the neutron EDM 
experiment. 

In the SU(5) SUSY GUT with right-handed neutrinos, the new flavor mixing in the 
right-handed squark sector is induced by the MNS matrix and GUT interactions. SUSY 
contributions to Cjl and Cjr can be as large as Cjj^. The EDM constraints are also strong. 
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FIG. 4: Wilson coefficients (a) C-jl and (b) C^r normalized by the SM value of C-jl- 



21 





mSUGRA 


SU(5) SUSY GUT 


U(2) 


degenerate 


non-degenerate 


C7,8L 


- 


- 








- 












- 








- 



















TABLE IV: Possible SUSY contributions to Wilson coefficients C7's and Cg's and Mi2(-Bs) in each 
model. ^ means non-negligible deviation from the SM, and denotes large SUSY contributions. 

and the SUSY contribution to the phase of Cy^ cannot become large as in the mSUGRA. 
Since we have introduced no CP phase of the neutrino sector in this analysis, the SUSY 
contribution to the phase of C'j^ mainly comes from the GUT phases 0f . Note that 0f 's 
contribute only to off-diagonal elements of the right-handed down squark mass matrix, and 
thus, they do not affect the neutron EDM. The /i 67 constraint in the degenerate case is 
much stronger than that in the non-degenerate case as seen in Fig.^ Therefore, the allowed 
regions become much larger in the non-degenerate case. 

In the U(2) model, SUSY contributions to C7L and C7R can be large because of the 
existence of new flavor mixings in the squark sector. Though the new contribution to the 
phase of C7L is restricted by the neutron EDM constraint in this model, the restriction is 
weaker than that in the above two models because the phase of (mQ)23 is independent of 
the phase that contributes to the neutron EDM. Compared with the non-degenerate case of 
SU(5) SUSY GUT with right-handed neutrinos, |(m|,)23| is suppressed in the U(2) model, 
and the SUSY contribution to 1 6*7/2! is smaller. 

The SUSY contributions to CgL and Cg/j are similar to those to C7L and C7R in each 
model as seen in Fig. El because the flavor mixings and CP phases that determine the SUSY 
contributions to C^l and C^r are the same as those contribute to Csl and Cg/j- This means 
that the constraints on C7L and C-jb. from b ^ S'y branching ratio also restrict predicted 
values of Csl and Csr. In Table HVl we summarize possible SUSY contributions to Wilson 
coefficients C7's and Cg's and Mi2{Bs) in each models. 

In Fig. ini we plot A^p{B —>■ Xsj) and A^p{B Mg'-f) versus the gluino mass for 
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FIG. 5: Wilson coefficients (a) CgL and (b) Csr normalized by the SM value of CgL- 
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tan/3 = 30. A'^'p{B X,7) essentially comes from the imaginary part of the interference 
terms between Cjl and C2L, and Csl and C2L, because C2R is negligible. Therefore, SUSY 
contributions to A^plB Xs'j) is constrained by the neutron EDM, and |yl^p(i? -^s7)| 
is at most ~ 1% in the mSUGRA and the SU(5) SUSY GUT with right-handed neutrinos. 
A^p{B Xs'y) \ can be as large as 3% in the U(2) model. The SM prediction is about 0.5% 
Iq . On the other hand, the mixing induced CP asymmetry in Ms'j depends on Cjl 

and Cjji. Although the SUSY contribution to Cjl cannot be negligible in the models that we 
are studying, the deviation from the SM essentially comes from Cjji. Therefore, the SUSY 
effect can become larger in the SU(5) SUSY GUT with right-handed neutrinos and in the 
U(2) model compared with the mSUGRA model. In the mSUGRA model, M^7)| 
is at a level of 1%, which is similar to the value of the SM In the non-degenerate case 
of the SU(5) SUSY GUT with right-handed neutrinos, {A'^plBfj Ms'~f)\ can be maximal, 
while in the degenerate case, {A'^p^Bd Mg'y)] can be as large as 0.1. In the U(2) model, 
we find that \Aqp{B Ms'y)\ could be as large as 0.5. 

In Fig. m we show the correlation between y4™p(i?d <t>Ks) and A^^p[Bd J/ipKs) 
for tan/? = 30. In the SM, ^ ^Ks) = Acp(^rf J/'ipKs) is satisfied. As 

mentioned in Sec. IIIIl the SM contribution is dominant in Cll, C^lr^ and C^^p due to the 
QCD correction between the electroweak scale and the bottom mass scale. Thus, SUSY 
contributes to A'^p^Ba 4>Ks) mainly through Cg's. 

In the mSUGRA, we see that the SM relation A^'p{Bd (f)Ks) = A^'f{Bd J/^Ks) 
approximately holds, and the deviation from the SM in A^^p{B J/ipKs) is less than 10% 
as seen in Ref. 0]. Therefore, y4™p(i? 4>Ks) is almost the same as that in the SM. 

In the non-degenerate case of the SU(5) SUSY GUT with right-handed neutrinos, 
A^^p{B(i —>■ (pKs) can substantially differ from the value in the SM and may be smaller than 
0.1, because of the large SUSY contributions to Csr- On the other hand, in the degenerate 
case, the SM relation A'^'piBd (pKs) = Acp(^d ^ J/^^Ks) is satisfied. Accordingly, the 
value of A^^p{Bd 4>Ks) is restricted by the experimental result on A^p{Bd J/ipKs). 

In the U(2) model, A^p{B (pKs) can deviate from the SM prediction because of SUSY 
contributions to Csl and Csr. The experimental result of A^p{Bd J/ipKs) implies that 
^ (t)Ks) lies between 0.3 and 1.0. 

In Fig. ISl we show y4™p (5^ <pKs) as a function of the gluino mass. In the mSUGRA 
and the degenerate case of the SU(5) SUSY GUT with right-handed neutrinos, the SUSY 
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FIG. 7: The correlation between the mixing-induced CP asymmetries in Bj 4>Ks and 
J/ipKs- The vertical and horizontal dotted lines show the la ranges of experimental values. In 
this plot, the experimental constraint of A™p(Sd — J/il)Ks) is not imposed. 

effect is almost negligible, and we see virtually no dependence on the gluino mass. On the 
other hand, in the degenerate case of the SU(5) SUSY GUT with right-handed neutrinos 
and the U(2) model, the SUSY contribution is significant, in particular for the smaller mass 
of the gluino. 

Recently, it has been pointed out that the contribution of the chromo-EDM of the strange 
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FIG. 8: The mixing-induced CP asymmetry in B^i — > (pKs as a function of the gluino mass. 

quark to the EDM of ^^^Hg is devastating in SUSY models with a large 2-3 mixing in the 
right-handed squark sector, and that a large deviation of A^p{Bd — » (t>Ks) from the SM 



prediction is unlikely provided that the experimental upper bound 5l| on the EDM of 
'»»Hg . imposed Q. Although the theoretical estimate M of the '»Hg EDM due to the 

chromo-EDM of the light quarks suffers from relatively large uncertainties, it is probable 
that the ^^^Hg EDM constraint severely restricts flavor mixings and CP violation in some 
SUSY models, and thus constricts their flavor and CP signals such as A^p{Bd (pKs) in 
general. 
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TABLE V: Significance of SUSY contributions to the CP asymmetries in eacli model, means 
non- negligible deviation from the SM, and ^J^J means large SUSY contributions. 

Among the models considered in the present work, effects of the ^^^Hg EDM constraint 
are non-neghgible in the non-degenerate case of the SU(5) SUSY GUT with right-handed 
neutrinos and the U(2) model. We show the correlation between the ^^^Hg EDM and 
A^^p{B(i (pKs) in these models in Fig. IHl In our numerical calculation of the ^^^Hg 
EDM, contributions of all the three light flavors are included. Fig. El illustrates how a flavor 
signal is tightened if the ^^^Hg EDM constraint is applied. In Fig.^[ the correlation between 
A'^^{Bd J/ipKs) and A'^'^{Bd (f)Ks) under the constraint of the ^^^Hg EDM is shown. 
Comparing this figure with FigEl we see that the ^^^Hg EDM constraint certainly restricts 
the possible deviation from the SM prediction. At the same time, however, our detailed 
calculation shows that there still remains some parameter region in which A'^p^Bd (pKs) 
is significantly different from A^'^{Bfi —>■ J/tpKs) in these models. A similar argument is 
applied to both A'^'p{B X,7) and A^^{B M^-f). 

In Table El we summarize the significance of SUSY contributions on the CP asymmetries 
that we have considered. In this table, we see the possibility to distinguish the three models 
in B experiments. 

As we stressed in the introduction, the purpose of this work is to demonstrate that 
identifying patterns of deviations from the SM predictions is useful to distinguish different 
origins of the SUSY breaking sector. From this point of view, combined with the analysis 
of the unitarity triangle|4|, we can make the following observations: 

• Deviations from the SM predictions in the unitarity triangle and rare decays are small 
in the mSUGRA model, except for some sizable contributions in the direct CP violation 
in the 6 — 57 process. Note that this conclusion may not hold in a particularly large 
value of tan /3 ~ 60 due to the Higgs exchange effects . 
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• The pattern of the deviations from the SM depends on the right-handed neutrino 
mass matrix in the SU(5) SUSY GUT with right-handed neutrinos. In the degenerate 
case, flavor mixing signals between the 1-2 generations become large. This appears as 
inconsistency between the measured value of ex and the B meson unitarity triangle, 
although the unitarity triangle is closed among B meson observables. The rare decay 
processes induced by the h-s transition do not show large deviations, but the branching 
ratio of /i — > 67 process can be just below the present experimental bound. This is 
expected to be a generic feature of SU(5) SUSY GUT with right-handed neutrinos. 

• In a specific parameter choice of the "non-degenerate" case, in which the — > 67 
constraint is relaxed, the flavor signals between 2-3 generations are expected to be 
sizable. This includes the mixing-induced CP asymmetry in B^ —>■ M^j and B^ — »• 
(pKs- The direct CP asymmetry in the b ^ s'y process, on the other hand, does not 
show a large deviation. 

• Various new physics signals in the consistency test of the unitarity triangle and rare 
decay process are expected in the MSSM with U(2) flavor symmetry. 

In this way, we can expect different sizes and patterns of new physics signals in the above 
models. These are crucial in pointing toward a specific model from flavor physics. 

V. CONCLUSIONS 

In order to seek the possibility to distinguish different SUSY models with B physics 
experiments, we have studied rare B decays related to the b s transition combining with 
the unitarity triangle analysis in three SUSY models. These models, namely the mSUGRA, 
the SU(5) SUSY GUT with right-handed neutrinos, and the U(2) flavor symmetry model, 
are different in character with respect to flavor structures of their SUSY breaking sectors. 
We have considered two different cases in regard to the mass spectrum of the right-handed 
neutrinos in the SU(5) SUSY GUT with right-handed neutrinos. 

In the unitarity triangle analysis, we have studied consequences of SUSY to A^p^B^ 
J/ipKs), AitlbJ ArriB^, and 03 . Our results are summarized in Table ITTl and Fig. [21 It could 
be possible to distinguish the three models by precisely measuring AniB^ and 03 in future 
B experiments. 
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As for rare B decays, we have explored SUSY effects to the direct CP asymmetry in 6 — >■ 
s'-f, the mixing induced CP asymmetry in Mg'-f, and the CP asymmetry in B^^ <pKs 

in the three models. The results are summarized in Tables IIVI and |^ Table IIVI shows the 
relative importance of SUSY contributions to the theoretically interesting Wilson coefficients 
related to the 6 — > s transitions and the Bg-Bg mixing amplitude. The significance of SUSY 
effects to the CP asymmetries is indicated in Table IVl 

The new flavor signals in the mSUGRA and the degenerate case of the SU(5) SUSY GUT 
are relatively limited in the 6 — > s rare decays considered in the present work. To detect 
these signals, typically a few percent, we may need an ultimate B experiment. 

On the other hand, the non-degenerate case of the SU(5) SUSY GUT exhibits quite 
attractive flavor signals in Bd M<j7 and Bd 4>Ks as seen in Table |3 We have also 
observed that the U(2) model predicts significant deviations from the SM in the 6 — > s rare 
decays as well as the unitarity triangle analysis. So far, both Belle and BaBar experiments 
have collected copious B decays, and they are expected to go well continuously. Thus, more 
Bd —>■ (pKs and related events will be obtained in near future. Moreover, both KEK and 
SLAG plan to upgrade their B factories. Therefore, the above flavor signals may well be in 
the reach of the present and foreseeable future B experiments. 

Combining the above observation with the results in our previous work, we conclude that 
the study of the unitarity triangle and rare B decays could discriminate several SUSY models 
that have different flavor structures in their SUSY breaking sectors. Such a study will play 
important roles, even if SUSY particles are found at future experiments at the energy frontier 
such as LHC. Although the spectrum of SUSY particles will be determined at LHC and a 
future e'^e~ linear collider, most of information concerning the flavor mixing of the squark 
sector is expected to come from the super B factory and hadron B experiments. Since the 
flavor structure of the SUSY breaking provides us with an important clue to the origin of 
the SUSY breaking mechanism and interactions at very high energy scales, B physics will 
be essential for clarifying a whole picture of the SUSY model. 
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